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We study the performance of initial product states of n-body systems in generalized quantum metrology 
protocols that involve estimating an unknown coupling constant in a nonlinear fc-body (k <C n) Hamiltonian. 
We obtain the theoretical lower bound on the uncertainty in the estimate of the parameter. For arbitrary initial 
states, the lower bound scales as l/n k , and for initial product states, it scales as l/n fe-1 / 2 . We show that 
the latter scaling can be achieved using simple, separable measurements. We analyze in detail the case of a 
quadratic Hamiltonian (k = 2), implementable with Bose-Einstein condensates. We formulate a simple model, 
based on the evolution of angular-momentum coherent states, which explains the 0(n~ 3 ^ 2 ) scaling for k — 2; 
the model shows that the entanglement generated by the quadratic Hamiltonian does not play a role in the 
enhanced sensitivity scaling. We show that phase decoherence does not affect the 0{n~ i ^ 2 ) sensitivity scaling 
for initial product states. 



I. INTRODUCTION 

Parameter estimation is a fundamental physical task. It typ- 
ically involves picking a physical system whose state, through 
its evolution, depends on the value of the parameter. In most 
quantum metrology schemes [1-18], this system, which we 
call the "probe," is a composite made up of n elementary 
quantum constituents. The influence of the unknown parame- 
ter 7 on the probe is described by an n-body Hamiltonian 



i? 7 = -fH 



(1.1) 



in which 7 appears as a coupling constant and H is a dimen- 
sionless coupling Hamiltonian (we use units with % = 1, so 
7 has units of frequency). The precision with which 7 can 
be determined depends on the initial state of the probe, the 
nature of the parameter-dependent Hamiltonian, and the mea- 
surements that are performed on the probe to extract informa- 
tion about the parameter. Other factors, such as decoherence 
in the probe [3, 19-21], also have an effect on the achievable 
sensitivity. 

The appropriate measure of the precision with which 7 can 
be determined is the units-corrected mean-square deviation of 
the estimate 7 0St from the true value 7 [22, 23]: 
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This estimator uncertainty is inversely proportional to the dis- 
placement in Hilbert space of the state of the probe corre- 
sponding to small changes in 7. The fundamental limit on 
the precision of parameter estimation, 



£7 > 
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(1-3) 
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called the quantum Cramer-Rao bound (QCRB) [22-25], is 
an expression of the maximum amount the state can change 
under the evolution due to H 1 . In Eq. (1.3), v is the number 
of trials with independent, identical probes, t is the time for 
which each probe evolves under H 1 , and AH = ((H 2 ) — 
(H) 2 ) 1 / 2 denotes the uncertainty in H for each probe, which 
does not change under the evolution due to iJ 7 . The QCRB 
is independent of the choice of estimator and is achievable 
asymptotically in the limit of a large number of trials, pro- 
vided the initial state of the probe is a pure state. If the initial 
state of the probe is not pure or if nonunitary processes de- 
stroy the purity of the initial state, the bound (1.3) is not tight, 
and a stricter version of the QCRB, given in [22-25], can be 
used, but we have no need for this stricter bound in this paper. 
The uncertainty in H is bounded above by 



AH < 



(1.4) 



where A max and A m i n are the maximum and minimum eigen- 
values of H. The difference between the largest and least 
eigenvalues, denoted by 



ii? II = A max — A n 



(1.5) 



is an operator semi-norm of H . Using this semi-norm, we can 
write a state-independent version of the QCRB [22] : 
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This bound can be achieved by using the initial state 
(|A max ) + I A m i n )) / %/2, which evolves after time t to 
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Measurement in a basis that includes the states |±) = 
(I A max ) ± lAmin))/-^ yields outcomes ±1 with probabilities 
p+ = cos 2 (||i7|| 7 i/2) andp_ = sin 2 ( 1 1 1 1 ^£/2) . Letting a 
denote an observable with the outcome values ±1, we have 
(a) = cos(||£f|| 7 i) and Act = |sin(||i?||7f)|. An appropriate 
estimator is defined in terms of the mean of the outcomes, 



1 " 

-^cr fe = cos(||i?||7e St i) 



k=l 

= oos(||ff||7t) - \\H\\ tsin(||H|| 7 t)(7 ert - 7) , 

(1.8) 

where the second expression is the linear approximation to the 
relation between 7 est and the mean of the outcomes, holding 
statistically in the limit of a large number of trials, specifically, 

v ^> tan 2 (||iJ||7i). Now it is easy to see that (7 es t) = 7 an d 
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showing that the bound (1.6) can be achieved and thus making 
it the fundamental limit to quantum metrology. 

The ^/v factor in Eqs. (1.3) and (1.6) is the well-understood 
statistical improvement available from averaging over many 
probes. In the remainder of the paper, we sometimes do not 
include this factor explicitly, referring to the remaining term 
on right-hand side of Eq. (1.3) or Eq. (1.6) as the QCRB, al- 
ways remembering, of course, that generally the bound can 
only be achieved asymptotically in the limit of large v. 

It is clear from Eqs. (1.3) and (1.6) that strategies for im- 
proving the precision in estimating a parameter include chang- 
ing the initial state of the probe, the coupling Hamiltonian H, 
or both. The case that has received most attention in the past 
is the one in which the probe constituents are coupled inde- 
pendently to the parameter: 



(1.10) 



Here hj is a single-body operator acting on the jth constituent 
of the probe (hence, all these operators commute). In this 
case, the QCRB (1.6) scales like <5 7 = 0(l/n), a scaling 
known as the Heisenberg limit [1, 10]. This scaling outper- 
forms that attainable with classical statistics, which goes as 
5 j = 0(1/ y/n), a scaling known as the standard quantum 
limit or, sometimes, as the shot-noise limit. The 1 / ' yfn scaling 
is the optimum sensitivity allowed by the QCRB (1.3) when 
the coupling Hamiltonian has the form (1.10) and the opti- 
mization of the probe state is restricted to product states. 

Achieving the quantum-enhanced Heisenberg sensitivity 
with the linear coupling Hamiltonian (1.10) requires the probe 
to be initialized in a highly entangled state, which is a 
formidable challenge using current technology [1]. Much 



progress has been made in preparing such states toward pre- 
cisely this end, e.g., by using measurement-induced squeez- 
ing [26], but it is still currently infeasible to use Heisenberg - 
limited experiments to outperform the best measurements op- 
erating at the standard quantum limit, although Heisenberg- 
like scalings have been achieved in related serial protocols 
that involve repeated interactions with a single constituent 
and thus do not involve entanglement [27]. Practical propos- 
als [28] for reaching the Heisenberg limit have also been made 
in the context of measurements on a harmonic oscillator pre- 
pared in a state that displays sub-Planck phase-space struc- 
ture [29]. 

More general families of Hamiltonians, containing nonlin- 
ear couplings of the constituents to the parameter, in contrast 
to the independent, linear coupling of Eq. (1.10), can perform 
better than the 1 /n scaling of the Heisenberg limit, while re- 
specting the QCRB [6, 9, 11, 12, 14, 15, 18]. In particular, 
when the coupling Hamiltonian has symmetric A: -body terms 
in it, it is possible to achieve the scaling Sj = (D(n~ k ), as 
shown in [12]. This (D(n~ k ) scaling requires entangled input 
states, but we will show here that the optimal scaling with ini- 
tial product states is 0(n~ k+1 / 2 ). Thus scalings better than 
the 1 jn Heisenberg scaling are possible for k > 2, even with 
initial product states for the probe, a result found for k = 2 
in [6, 9, 18]. 

In this paper we investigate the theoretical and practical 
bounds on precision in the generalized quantum metrology 
scheme introduced in [12], which allows for nonlinear cou- 
plings of the probe constituents to the parameter. In particular, 
we study the nonlinear coupling Hamiltonian [12] 
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For simplicity, we assume that the probe constituents are iden- 
tical and that the single-body operators hj are the same for all 
the constituents. We review in Sec. II the optimal precision 
that can be achieved when the probe can be prepared in any 
initial state, particularly, entangled states of the probe con- 
stituents, but our emphasis in this paper is on the precision 
that can be attained when the initial state is restricted to be a 
product state. 

In Sec. Ill we show that the optimal precision with product- 
state inputs scales as (D(n~ k+1 / 2 ), and we find the corre- 
sponding optimal input state for the probe constituents. The 
sensitivity that can be achieved in practice depends, as men- 
tioned above, on the measurements that are performed on the 
probe to extract information about 7. By analyzing the short- 
time limit, we show in Sec. IV that when the probe is ini- 
tialized in a product state, simple, separable measurements on 
the probe constituents can achieve the optimal sensitivity. The 
conclusion, reached in [6, 9, 18] for specific cases, is that scal- 
ings better than the standard quantum limit — indeed, better 
than the Heisenberg limit — can be had without the need to in- 
vest in the generation of fragile initial entangled states. Our 
scheme thus circumvents a major bottleneck in attaining, in 
practice, scalings superior to 1/ y/n. 

In Sec. V, we analyze in detail a quadratic coupling Hamil- 
tonian (k = 2) for effective qubits, with hj — Zj/2, where Zj 
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is the Pauli z operator for the jth qubit. This case can be im- 
plemented in Bose-Einstein condensates [15, 18], as was sug- 
gested in [12], and in fermionic atoms in optical lattices [15]. 
We show that the optimal sensitivity for input product states, 
scaling as 0(n~ 3 / 2 ), can be achieved by using the optimal 
product-state input found in Sec. Ill and making separable 
measurements of equatorial components of the total angular 
momentum of the effective qubits. Moreover, we show that, 
as was found independently by [18], the 0(n~ 3 / 2 ) scaling is 
attainable with these measurements starting with almost any 
state of the constituents except the equatorial states, which 
were the subject of the analysis in [15]. We formulate a sim- 
ple model, based on the evolution of angular-momentum co- 
herent states, that explains the origin of the 0(n~ 3 / 2 ) scaling. 
The model indicates that the entanglement generated by the 
quadratic Hamiltonian does not play a role in the enhanced 
sensitivity, and it suggests that, unlike protocols based on the 
use of entangled inputs, the product-state scheme should not 
be extremely sensitive to decoherence. We verify this sugges- 
tion by a brief analysis of the effect of phase decoherence in 
Sec. VB. 

Section VI concludes with a brief summary of our results, 
including an extension of the coherent-state model to arbitrary 
fc, and a discussion of our perception of the field of quantum- 
enhanced metrology. 



II. GENERALIZED QUANTUM METROLOGY AND THE 
QUANTUM CRAMER-RAO BOUND 

Attaining the QCRB (1.6) requires using an appropriate ini- 
tial state and making appropriate measurements to extract the 
information about 7, but changing the optimal scaling with 
n requires changing the dependence of the coupling Hamil- 
tonian H on n. This can be done by replacing the linear 
coupling Hamiltonian (1.10), which has just n terms, with 
the nonlinear Hamiltonian (1.11). The coupling Hamilto- 
nian (1.11) describes a system with symmetric fc-body cou- 
plings, including self-interactions, and it has n k terms. For 
example, if the constituents are spin-i particles and the oper- 
ator hj is the z-component of the the jth particle's spin, then 
H describes a coupling of the parameter to the fcth power of 
the z-component of the total angular momentum. 

The eigenvectors of H are products of the eigenvectors of 
the hj&. The eigenvectors can be labeled by a vector of single- 
body eigenvalues, I A) = |Ai,...,A„). The corresponding 
eigenvalues of h are given by the polynomial 



U k (\) = (X\H\\)= HTa, = ]T A 



(2.1) 

which is symmetric under permutation of its arguments and 
is known in the mathematical literature as the fcth-degree ele- 
mentary symmetric polynomial (on n variables). 

To calculate the QCRB (1.6) for the fc-body coupling 
Hamiltonian (1.11), we have to calculate the maximum and 
minimum eigenvalues of H in terms of the eigenvalues of 
the single-body operators hj, the total number of constituents 



n, and the degree of the coupling fc. Let A E 
the largest and smallest eigenvalues of hj. We consider four 
cases. 

1. fc odd. The largest (smallest) eigenvalue of H is 

A max = (nA max ) fc [A min = (nA min ) fe ], correspond- 
ing to the eigenvector |A max ) = |A max , . . . , A max ) 

(|A m in) — |A m i n , . . . , A m i n )). 



2. k even, A„ 
Case 1. 



> 0. The same conclusions apply as in 



3. fc even, A max < 0. The same conclusions apply as 
in Case 2, except that the roles of A max and A m ; n are 
reversed: the largest (smallest) eigenvalue, A max = 
(nA min ) fc [A min = (nA max ) fe ], corresponds to the 
eigenvector that has every constituent in the state | A m i n ) 
(|A max )). 

4. k even, A m i n < < A max . Since all the eigenvalues of 
H are nonnegative, the maximum eigenvalue is A max = 
HA| max ) fe , where |A| max = max{|A max |, |A mi „|}, 
corresponding to all the constituents being in either 
|A max ) or I Amin)- The minimum eigenvalue comes 
from the string A, perhaps containing all eigenvalues, 
that makes (A) as close to zero as possible. 

In Cases 1-3, the QCRB (1.6) takes the form 

1 



h > 



^ fe |A maX - A mi 



(2.2) 



displaying the 0(n~ k ) scaling found in [12]. 

Case 4 requires further discussion regarding A m in- We can 
bound Amin from above by considering strings A that contain 
only A max and A m i n . If A max appears a fraction p of the time, 
the corresponding eigenvalue is 



n fe (A) = [npA max + n(l -p)A min ] fc ■ 

This eigenvalue can be minimized by making pA max 
p)A m in as close to zero as possible, i.e., by choosing 



np - 



n|A n 



(2.3) 
(1- 

(2.4) 



where \x\ denotes the nearest integer to x and \\h\\ = A max — 
Amin is the semi-norm of the single-particle operators. The 
resulting eigenvalue can be written as Hfe(A) = ((5||/i||) fe , 
where <5 < 1/2 is the magnitude of the difference between 
«|A m i n |/||/i|| and the integer closest to it. The minimum 
eigenvalue can be written as 



A min = [(<y-e)||/i||: 



(2.5) 



where e (0 < e < 1/2) accounts for the fact that strings 
containing other eigenvalues of the hjS can generally make 
A m in smaller. If the constituents are qubits, e = 0. Since 

IHI/|A| max < 2, we have 



A n 



A n 



< n 



(2.6) 
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The upshot is that for even k and A m j n < < A max , the 
QCRB (1.6) is given by 



' tn k \X\ k 1 - A • /A ' 

I I max n mm/ n max 

Thus the symmetric fc-body coupling leads to a QCRB scaling 
as 0(n~ k ) in all four cases. 

A closely related fc-body coupling Hamiltonian is the same 
as Eq. (1.11), except that the self-interaction terms are omit- 
ted, which might be more appropriate in some physical situ- 
ations. We analyze this alternative fc-body coupling Hamilto- 
nian in App. A and show that it also leads to a 0(n~ k ) QCRB 
scaling when arbitrary input states are allowed. 

Luis and collaborators [6, 9] showed that 0(n~ 2 ) scalings 
can be achieved in principle using a Kerr-type optical non- 
linearity, and Luis generalized these results to optical non- 
linearities of arbitrary order in [14]. Reference [17] pro- 
poses a method for synthesizing a quadratic (k = 2) Hamilto- 
nian from a linear Hamiltonian by passing a light beam twice 
through an atomic medium and finds a <D(n~ 2 ) scaling for this 
method. 



fewer terms. The first sum in Eq. (3.2), in which the terms 
have no duplicate factors, has nl/(n — k)\ = 0(n k ) terms, 
and the second sum, in which one factor is duplicated in each 
term, has n\/(n — k — 1)! = ©(n*^ 1 ) terms. The binomial 
coefficient multiplying the second sum accounts for the num- 
ber of ways of choosing the factor that is duplicated. The next 
sums in the expansion, involving terms with factors and 
h 2 hf, have n\/{n — k — 2)! = (D(n k ~ 2 ) terms. These expan- 
sions require that n > k, which we assume henceforth, and 
the scalings we identify further require that n^> k. 

Given the expansion (3.2), the expectation value of H has 
the form 

(H)= E (h ai ) ■ ■ ■ (h ak ) 

(a 1 ,...,a k ) 




E (ha 1 )---{ha k _ 2 )(hl k _ 1 ) 



+ 0(n k - 2 ). (3.3) 
The expression for {H 2 ) follows by replacing k with 2k: 



III. ATTAINABLE PRECISION WITH PURE PRODUCT 

STATES 

A. General bound 

The QCRB (1.6) gives the best possible measurement pre- 
cision, but can only be achieved for an optimal probe initial 
state, i.e., one of the form (|A max ) + e 1 ^ |A m ; n ))/\/2, which 
the results of Sec. II show is typically highly entangled. In this 
section we obtain lower bounds on £7 in the situation where 
the initial state is a pure product state, 

|*o) = hM®---®hM ; (3.1) 

for this purpose, we start from the state-dependent 
QCRB (1.3). Since all the one-body operators hj in the cou- 
pling Hamiltonian are assumed to be identical, it is reason- 
able to expect that the optimal initial product state will have 
all constituents in the same state, but we do not assume this 
at the outset, instead allowing its moral equivalent to emerge 
from the analysis. 

The trick to evaluating AH is to partition the unrestricted 
sum in Eq. (1.1 1), in which terms in the sum contain different 
numbers of duplicate factors, into sums such that each term 
has the same sort of duplicate factors. Thus we write 

H = E h ai - ■ - h ak 

(ai,...,a k ) 




(3.2) 



where a summing range with parentheses, (ai, . . . , a;), de- 
notes a sum over all /-tuples that have no two elements equal. 
The two sums in Eq. (3.2) are the leading- and subleading- 
order terms in an expansion in which successive sums have 



(H 2 )= E (h ai )---(h a2k ) 

(a 1 ,...,a 2k ) 

+ (?) E (h ai )---(h a2k _ 2 )(hl 2k _ 1 ) 

(ai ,...,a 2 fc-i) 

+ 0(n 2k - 2 ). (3.4) 

The rest of the analysis is based on an artful switching be- 
tween restricted and unrestricted sums. By changing the ini- 
tial sum in Eq. (3.3) to an unrestricted sum, we can rewrite 
(H) to the required order as 

(m= E CO ■■■<&.*> 

ai,...,a k 




E (^ 1 >"-<^_ 2 )A^_ 1 



+ 0(n k - 2 ). (3.5) 

Squaring this expression and changing the unrestricted sums 
back to restricted ones, again keeping only the leading- and 
subleading-order terms, gives 

(Hf= E (h ai )---(h a2k ) 

(ai,...,a 2k ) 




+ 2Q) E (ha 1 )---(K 2k _ 2 )Ah 2 ak _ 1 

(ai,...,a 2 je_i) 

+ 0(n 2k - 2 ). (3.6) 
We can now find (AH) 2 by subtracting Eq. (3.6) from 
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Eq. (3.4): 

(AH) 2 = k 2 £ (h ai )--(h aah _ a )Ahl k _ 1 

(a 1 ,...,a 2k -i) 

+ 0(n 2k - 2 ) 

/ ™ \ 2(fc-l) , n \ 

(3.7) 

In the final form, we take advantage of the fact that in the now 
leading-order sum, we can convert the restricted sum to an 
unrestricted one. 

To make the QCRB (1.3) as small as possible, we need to 
maximize the variance (AH) 2 of Eq. (3.7). We can imme- 
diately see that for fixed expectation values (hj), we should 
maximize the variances Ah 2 , and this is done by using for 
each constituent a state that lies in the subspace spanned by 
|A max ) and |A m i n ). Letting pj be the probability associated 
with | A max ) for the jth constituent, we have 

X 3 = (hj) = PjA max + (1 - Pj) Xmin = Amin + Pj\\h\\ , 

Ah 2 = PjA max + (1 - p j )X 2 nin - x 2 
= \\hfpj(l-pj) 

= (A max - Xj)(xj - A min ) . (3.8) 
Thus we should maximize 

/ n \ 2(fe— 1) n 
(AH) 2 = k 2 ( Xj ; j ^(A max - Xj)(Xj - A min ) 

(3.9) 

within the domain defined by A m j n < Xj < A max , j = 
l,...,n. 

Discarding potential extrema of (A_ff) 2 given by = 
J2j Xj, since these either are minima or lie outside the rel- 
evant domain, we find that the conditions for extrema of 
(AH) 2 imply immediately that Xj — x (and thus pj = p) 
for j = 1, . . . , n. Thus the optimal states in the initial product 
state (3.1) have the form 

\*l>j) = VP |Am aX ) + v/l - ~p | Amin) • (3.10) 

The only possible difference between the states for differ- 
ent constituents is in the relative phase between |A max ) and 

I A m in). 

Since the optimal constituent states live and evolve in a two- 
dimensional subspace, we can regard the constituents effec- 
tively as qubits, with standard basis states |0) = |A max ) and 
|1) = |A m i n ), serving as the basis for constructing Pauli oper- 
ators X, Y, and Z. Restricted to this subspace, the operator h 
takes the form 

h = A max |0> <0| + A min 1 1> < 1 1 = Al + \\h\\Z/2 , (3.11) 

where A = (A max + A m i n )/2 is the arithmetic mean of the 
largest and smallest eigenvalues of h. 



In the analyses in Sees. IV and V, we assume that all the 
constituents have zero relative phase (<pj = 0), giving an ini- 
tial state \V ) = \ipp)® n , where 

= e~^ Y ' 2 |0) = cos(/3/2) |0) + sin(/3/2) |1) , 
p = cos 2 (/3/2) = (l + cos/3)/2 . (3.12) 

Here we describe the one-body state in terms of the ro- 
tation angle [3 about the y axis that produces it from |0). The 
corresponding initial density operator is 

n 

Pfl = | * /3 ) | = (g) - (1 j + Xj sin + Zj cos (3) . (3 . 13) 

3 = 1 

The variance of H now takes the simple form 

(AH) 2 - fcV*- 1 (h) 2 ^ (Ah) 2 

= fcV*- 1 * 2 **" 1 )^ -x)(x- A min ) , (3.14) 

which leads, in the QCRB (1.3), to a sensitivity that scales as 
l/n fe_1 / 2 for input product states. This should be compared 
with the 0(n~ k ) scaling that can be obtained by using initial 
entangled states [12]. Notice that for k > 2, the 0(n^ k+1 / 2 ) 
scaling is better than the 1/n scaling of the Heisenberg limit, 
which is the best that can be achieved in the k = 1 case even 
with entangled initial states. 

The coupling Hamiltonian that has the self-interaction 
terms of Eq. (1.11) removed is analyzed in App. A. We show 
that for initial product states, this modified Hamiltonian has 
the same leading-order behavior in the variance of H; it thus 
has 0(n~ fc+1 / 2 ) scaling for initial product states and the same 
optimal product states as we now find for the coupling Hamil- 
tonian (1.11). 



B. Optimal product states 

The problem of finding the optimal input product state is 
now reduced to maximizing the 2fc-degree polynomial 

f(x) = x 2(fc - 1} (A max - x)(x - A min ) 

= a; 2(fe_1) (||/i|| 2 /4- (x - A) 2 ) (3.15) 

with respect to the single variable x = (h) on the domain 
Amin < x < A max . The condition for an extremum is 

= f'(x) = 2x 2k -*[(k-l){\\h\\ 2 /A-(x A) 2 ) 

- x(x - A)] . (3.16) 

We assume k > 2, because the k = 1 case is already well 
understood. For k = 1, there is a single maximum at x = 
A, corresponding to equal probabilities for |A max ) and |A m i n ) 
and to (AH) 2 =n||/i|| 2 /4. 

The polynomial / vanishes at x = A min and x = A max . We 
can make some general statements about the extrema of / in 
three cases. 
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1. If A m ; n < < A max , / has a minimum at x = 
and two maxima within the allowed domain, one at 
a positive x + > A and one at a negative x_ < A. 
The global maximum is at x + (x-) if |A| max = A max 

(|A| max — |A m i n |)- 

2. If A max > A m j n > 0, / has a maximum at x = 0, a 
minimum for a positive x_ < A m j n , and a maximum 
within the allowed domain at x + > A. Only the last of 
these lies in the relevant domain. 

3. If A m j n < A max < 0, / has a maximum at x = 0, a 
minimum for a negative x + > A max , and a maximum 
within the allowed domain at < A. Only the last of 
these lies in the relevant domain. 

These general observations are perhaps more enlightening 
than the form of the (nonzero) solutions of Eq. (3.16): 

"-('-s^s^ + O-i)'*"'- <3J7 > 

The ± here means the same thing as in the discussion of the 
three cases above. 

As k increases, x+ approaches A ma x, and x- approaches 
A m ; n . Indeed, as k — > oo, we have x + = (1 — 
l/2fc)A max , corresponding to p + = 1 - A max /2fc||/i|| 
and (AH) 2 = (k/2e)(n\ max ) 2k - 1 \\h\\, and x_ = (1 - 
l/2fc)A m ; n , corresponding to p_ = — A min /2fc||/i|| and 
(AH) 2 = (k/2e)(-n\^) 2k - 1 \\h\\. 

An important limiting case, not covered in the discussion 
above, occurs when A m j n = — A max . Then the maxima occur 
symmetrically at 

x± =±±\\h\\y/l -1/k, (3.18) 

corresponding to probabilities p± = ^ + x±/\\h\\ = |(1 ± 
\fl — 1/k) = 1 — p T and to sin j3± = \Jljk . The two 
maxima lead to the same variance, 

(AH) 2 = k(l - l/k) k - 1 n 2k - 1 (\\h\\/2) 2k , (3.19) 

thus yielding a QCRB 

2 fe-i i 
5l ~ k 1 /2(l-l/k)( k - 1 V 2 tn k - 1 / 2 \\h\\ k ' (3 ' 20) 

Of course, when A m ; n = — A max , we can always choose units 
such that A max = l/2(||/i| = 1), which means that the single- 
body operators are hj = Zj/2. It is this situation that we 
analyze in the remainder of this paper. 



IV. SEPARABLE MEASUREMENTS 

In the previous section we obtained the theoretical limits 
on the measurement uncertainty with symmetric fc-body cou- 
plings and initial product states for the probe. The theoretical 



bound is saturated by a measurement of the so-called sym- 
metric logarithmic derivative [22-25]; this measurement, in 
general, is entangled and depends on the value of the param- 
eter that we are attempting to estimate. In this section we 
show that for some Hamiltonians of interest, standard separa- 
ble measurements lead to uncertainties for small 7 that have 
the same scaling as the theoretical bounds. The restriction to 
small values of 7 is not a strong limitation, because we can 
always use feedback to operate in this regime, as we discuss 
in more detail in Sec. V A. 

We consider the special case in which the single-body op- 
erators are hj = Zj/2, leading to a coupling Hamiltonian 




(4.1) 



Here we introduce J z as the z component of a "total angular 
momentum" corresponding to the effective qubits. We assume 
an initial state of the form (3.13), and we let this state evolve 
for a very short time, i.e., <f> = yt -C 1. In the remainder of the 
paper, we often work in terms of the parameter <j> instead of 7. 
After the time evolution, we measure the separable observable 

Jy = J2 Y J- ( 4 - 2 ) 

3 

Over v trials, we estimate <j) as a scaled arithmetic mean of the 
results of the J y measurements. 

The expectation value of any observable at time t is given 

by 

(M) t = Tr (uH'IUpp) = (U^MU) , (4.3) 

where U = e~ lH -' 1 = e^ tH ^, and where we introduce the 
convention that an expectation value with no subscript is taken 
with respect to the initial state. For small <f), we have 

U^MU = M — i(f)\M, H] + 0(<t> 2 ) . (4.4) 

Thus the expectation value and variance of J y at time t take 
the form 

(Jy)t - ( J y) - i( t> (i J y> H }) + > ( 4 - 5a ) 

(Ajy) 2 = (AJy) 2 

-i4>((Jy (Jy))[Jy,H} + [Jy,H](jy (Jy))) 

+ 0(<f> 2 ) . (4.5b) 

The initial expectation value and variance of J y are those 
of an angular-momentum coherent state in the x-z plane: 

(J y ) = , (4.6a) 

(AJy)l = (j 2 y ) = \j2W) = ! l- ^ 6b ) 

3,1 

In evaluating the other expectation values in Eqs. (4.5), we 
can avail ourselves of the expansions used in Sec. Ill A, since 
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we are only interested in the leading-order behavior in n. To 
leading order, the coupling Hamiltonian has the form 



H =^k E Z ai ---Z ak +0(n k -i). 



(4.7) 



(01, •••>«&) 



In this section we use w to indicate equalities that are good to 
leading order in n. We can now write 

[ J ^]«WiE E [Yj,Z ai ---Z ak ] 



j — l (ai,...,Ofc) 



~9fcE E ' ' ' ^ai-i^a,Z ai+1 ■ ■ ■ Z a 



2 k 

ik 
2^ 



(=1 (ai,...,Ofc) 



E ^ a2 ' ' ' ^ ak ' 



(4.8) 



(ai,...,a fe ) 

from which it follows that 

E (X ai )(Z a2 )---(Z ak ) 

(a 1 ,...,a k ) 



(4.9) 



Elaborating this procedure one step further, we can show that 
to leading order in n, the expectation value in the second line 
of Eq. (4.5b) vanishes. Our results to this point are summa- 
rized by 



= 4>k(n/2) k sin f3 cos^ 1 /? + 0((f> 2 ) , (4.10a) 
(AJ„) t w + O(0 2 ) . (4.10b) 

If we let our estimator </> ost be the arithmetic mean of the v 
measurements of J y , scaled by the factor (d (J y ) t /dcb)^ 1 = 
l/k(n/2) k smj3cos k ~ 1 j3, we have 



(4>est) = 



(Jy)t 



d{J y ) t /d<t> 

1 (AJ„)« 
07|rf(J y ) t /#| 



+ 0W> 2 ), 



(4.11) 



y/vkn k 1 ' 2 sm/3| cos fc 

This scheme thus attains the 0{n~ k+1 / 2 ) scaling that is the 
best that can be achieved by initial product states. In an anal- 
ysis of optical nonlinearities of arbitrary order, Luis [14] re- 
ported finding this 0(n~ k+1 / 2 ) scaling. 

The minimum of 8<j), occurring when sin (3 = \j\jk, gives 
an optimal sensitivity 



^ M ^*Va(l -l/fc)(*-D/> „*-V> +°^' (4 ' 13) 

which is identical to the optimal QCRB sensitivity for initial 
product states. For k — 2, the case that is the subject of the 



next section, the two optimal values of j3 are (3 = ir/4 and 
(3 = 37r/4, and the sensitivity becomes 



1 



+ o{4>) . 



(4.14) 



Aside from showing that the QCRB scaling for initial prod- 
uct states can be achieved, the analysis in this section serves 
to illustrate how the product-state scheme works in a regime 
that has a singularly simple description. The J k coupling 
Hamiltonian induces a nonlinear rotation about the z axis, 
which rotates the probe through an angle ( J y ). j (J x ) w 



# {Jz) 



k-1 



This rotation induces a signal in J y of size 



0fc (Jx) (Jz) k ^ 1 , which is k (J z ) k ^ 1 times bigger than for 
k = 1, yet is to be detected against the same coherent-state 
uncertainty y / n/2 in J y as for k = 1. To take advantage 
of the nonlinear rotation, we can't make the J x lever arm of 
the rotation as large as possible, because the nonlinear rota- 
tion vanishes when the initial coherent state lies in the equa- 
torial plane. Nonetheless, we still win when we make the 
optimal compromise between the nonlinear rotation and the 
lever arm. The optimal compromise comes from maximizing 
(Jx) (Jz) k 1 7 which turns out to be exactly the same as find- 
ing the optimum in the QCRB analysis of Sec. Ill B because 
(X) = sin 13 = AZ. 

A more careful consideration of the terms neglected in this 
analysis suggests that, as formulated in this section, the small- 
time approximation requires that cj> <C l/n k ~ 1 . Nonethe- 
less, the analysis is consistent because <j> can be resolved more 
finely than this scale, i.e., 8<f> n k ~ 1 —0(\/^/n). This conclu- 
sion is confirmed by the more detailed analysis of the k = 2 
case in Sec. V. On the other hand, the simple model of 
coherent-state evolution, developed for k = 2 in Sec. V, sug- 
gests the description of the preceding paragraph can be ex- 
tended to much larger times. We return to this point in the 
Conclusion. 



V. SEPARABLE MEASUREMENTS FOR THE 
INTERACTION H y = 7 J z 2 

We focus now on the symmetric, k = 2 coupling Hamilto- 
nian 

^7=7^ = 7(^^/2) . (5.1) 



This is perhaps the most important example for practical 
applications of nonlinear Hamiltonians to quantum metrol- 
ogy [12], since it occurs naturally whenever the strength of 
two-body interactions is modulated by a parameter 7. As sug- 
gested in [12], one good place to look for this kind of cou- 
pling Hamiltonian is in Bose-Einstein condensates. Indeed, 
in [15, 18], it is shown how this Hamiltonian can be imple- 
mented using the internal atomic states of BECs. In analyz- 
ing the BEC scenario, Ref. [15] finds a sensitivity that scales 
as 0(1/ n) for separable measurements made on a probe 
that evolves from an an initial product state chosen to be an 
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angular-momentum coherent state in the equatorial plane. The 
results in the previous sections show that we should be able to 
improve this scaling to 0(n~ 3 / 2 ) through a wiser choice of 
the initial coherent state. In this section we analyze this situa- 
tion in some detail. 

We take the initial state of the n-qubit probe to be an 
angular-momentum coherent state that is at an angle from 
the z axis in the x-z plane. This state is obtained from the 
coherent state along the z axis, | J, 3) = |0}®", by a rotation 
through about the y axis: 



m=e-* J y\J,J) = (e 



if3Y/2 



|0» 



C3>>. 



(5.2) 



The rotation about the y axis and the nonlinear rotation under 
the interaction Hamiltonian (5.1) both leave the state in the 
(23 + 1) -dimensional subspace with angular momentum 3 = 
n/2, so we can use the basis | J, to) of J z eigenstates for this 
subspace, with m = — 3, . . . , 3. The initial probe state used 
in [15] is a special case, = —n/2. 

The state = X)m=-j dm \ J, m), can be expanded in 
the basis | J, to) using a reduced Wigner rotation matrix [30] 



d m = d J mJ {0) = {J,m\e~ i0Jv \J,J) 

[cos(/3/2)] J+m [sin(/?/2)] J - m 



(2J)! 



(J + m)!(J -m)V 



At time t = <t>/~f, the state of the probe becomes 



(5.3) 



\* fj (t)) = e-V-Z = ]T d rrfi-^ | J, m) . (5.4) 



m=—J 



A. Measurements 



We now look at the attainable measurement uncertainties 
using both J x = and J y — ^\ measurements 

on the final state of the probe. It turns out that 3 X and J y 
measurements are on nearly the same footing, with 3 y mea- 
surements being marginally better, for all except = n/2. 
For very short times, the superiority of J y measurements for 
7^ n/2 is clear from the analysis in Sec. IV, since the 
change in ( J y ) is linear in <f), whereas the change in ( J x ) is 
quadratic. What happens for longer times and for = n/2 
cannot be addressed by the short-time analysis in Sec. III. 
What we find in this section is that both J x and J y measure- 
ments can achieve the optimal scaling obtained in Sec. Ill . 
For — n/2, J y measurements provide no information about 
(j), but J x measurements achieve the 0(n^ r ) scaling found 
in [15]. These conclusions assume no decoherence, and in 
Sec. VB we explore the impact of decoherence on the abil- 
ity to achieve super-Heisenberg scalings with the symmetric, 
k = 2 coupling Hamiltonian. 

For measurements of J x or J y , the sensitivity is given by 



c,y — t ^ X 



(^J X ,y*)(p 
|rf(J XiJ/ ) /#| 



(5.5) 



(for this subsection, we revert to our practice of omitting the 
1 / ^/v statistical factor from our sensitivity formulas). 

The expressions needed to calculate 6(f) for J x and J y mea- 
surements are derived in Appendix B. These results are con- 
veniently expressed in terms of the raising and lowering oper- 
ators 



J± J x i iJy j 



(5.6) 



since we can write 



(J x ) rj> = Re((J + } 4> ) , (J y )^ = lm((J + ) 4> ) , (5.7) 

(Jlv), - \ (J+J- + J-J+U ± lM (4), ) - (5-8) 

where the upper sign in Eq. (5.8) applies to J x and the lower 
sign to J y . In Appendix B, we show that 



(J+)^ = Jsin/3(cos 4> + isin(j)coa 0) 2J 1 
= J S m0r 2J - 1 e^ 2J -^ 6 , 



3(23 - 1) 



-(3 + 3_+3_3 + ) 4> = 3+ 2 



sin 2 /? 



(5.9) 
(5.10) 



(Jl) 4> = J(2J 2 1} sin 2 /3(cos 20 + i sin 20 cos 0) 2 ^-^ 

= J(2J 2 ~ 1} sin 2 /?^- 7 - 1 ^- 7 - 1 ) , (5.11) 

where 

r = (1 - sin 2 0sin 2 /3) 1/2 , (5.12a) 

6 = tan" 1 (tan <j) cos 0) , (5 . 1 2b) 

and 

R = (1 - sin 2 20sin 2 /3) 1/2 , (5.13a) 

6 = tan~ 1 (tan20cos/3) . (5.13b) 

Plugging these results into Eqs. (5.7) and (5.8), we arrive at 

(J^ = Jsin^r^cospJ- 1)0] , (5.14) 

(3 y ) 4> = 3sm0r 2J - 1 sm[{23 - 1)9] , (5.15) 

and 

, j2 . 3 3(23- 1) . 2o 

x (l±i? 2 < J - 1 )cos[2(J-l)e]) . 

(5.16) 

In using these results in what follows, it is easier to deal di- 
rectly with the first forms in Eqs. (5.9) and (5.11) rather than 
working with the functions r, 9, R, and G. 

The expectation values (3 XtV ) ^ change sign when <j> ad- 
vances by n. This means that their squares and absolute val- 
ues, which are all that appear in the sensitivity (5.5), are peri- 
odic with period n. The second moments ( Jj y ) are periodic 

with period n/2. The upshot is that the uncertainties A3 XtV 
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and the precision 8<ft x . y are periodic with period it. This ir- 
periodicity is a consequence of periodic revivals in the evolved 
state | */*(*)). 

The main features of the sensitivity 5(f) for measurements of 
J x and J y can be gleaned from Fig. 1. It is clear from these 
plots that the best sensitivity is achieved when (ft is near zero 
and also, because of the periodicity of 5(f), when (ft is near qir, 
for q any integer. 

When J is large, we can develop a good approximation for 
the entire region of high sensitivity, where (ft is small, by writ- 
ing 

(cos0 + *sin</> cos/3) 27 - 1 ~ e 2lJ * cosf) e" 14 ' 2 sin ' , 

(5.17a) 

(COs20 + isin2 ( />COS/3) 2 ( J - 1 ) ~ cos fi^ij^, 2 si„ 2 /3 _ 

(5.17b) 

These approximations are good to second order in (ft in the 
exponent. When (ft is near qir, the same approximations can 
be had by replacing (ft with (ft — qir. The complex exponentials 
give rise to rapidly oscillating fringes in {J x . y ) ^ and (J x<y ) ^, 
with periods ~ 1/ J cos (3; the slower Gaussian envelopes take 
these expressions to zero when (ft is a few times | sin. f3\/yfj. 

It is not hard to work out the sensitivity in this approxi- 
mation, but the formulas are sufficiently messy that they are 
little more illuminating than the exact expressions. We can, 
however, develop a very simple, yet instructive picture of the 
fringes by keeping them, but assuming that (ft is small enough 
that the Gaussian envelopes have yet to become effective, i.e., 
\/~J (ft sin (3 is somewhat smaller than 1 . In this approximation, 
the fringes are uniform in (ft, and we obtain 



where q and s are integers. At these operating points, the 
sensitivity for both measurements becomes 



(Jx)^ — Jsin/3cos(2J^cos/3) , 
{Jy), — Jsin/3sin(2J(/>cos/3) , 



(5.18) 
(5.19) 



(AJ,) 2 ~ -[1 - sin 2 /?cos 2 (2J</>cos/?)] , (5.20) 



(AJ,) 2 ~ l[l - sin 2 /3sin 2 (2J0cos/3)] 

These lead to sensitivities 

2 1 1 — sin 2 /3cos 2 (2J0cos/3) 
x ~ 2J3 sin 2 2/3sin 2 (2J0cos/3) ' 
2 _ 1 1 - sin 2 (3 sin 2 {2 J (ft cos /3) 
Yy ~ 2J 3 sin 2 2/3 cos 2 (2J0 cos (3) 



(5.21) 

(5.22) 
(5.23) 



Within this uniform-fringe approximation, the best sensi- 
tivities are achieved at the troughs of the fringes: the best op- 
erating points are, for J x measurements, 



qir 



(s + 1/2)tt 
2 J cos (3 



and for J y measurements, 



= qn + 



2 J cos [3 ' 



(5.24) 



(5.25) 



"x,y — 



V2J 3 / 2 |sin2/3| ' 



(5.26) 



which takes on its optimal value, l/^/2J z / 2 = 2/n 3 / 2 , when 
(3 = 7r/4 or (3 — 37r/4, in agreement with the analyses in 
Sees. Ill and IV. This 0{ J -3 / 2 ) scaling has been found inde- 
pendently by Choi and Sundaram [18] and for nonlinear opti- 
cal systems by Luis and collaborators [6, 9]. For (3 = 7r/4 and 
J = 2 500, Fig. 2 plots the central fringes of the approximate 
sensitivities (5.22) and (5.23) and compares them with the ex- 
act sensitivities and the Gaussian approximation for measure- 
ments of J x and J y . 

As we show in App. B, within the uniform-fringe approx- 
imation, the evolved state (5.4) is an angular- momentum co- 
herent state that makes an angle (3 with the z axis and that 
rotates around the z axis with angular velocity 27 J cos (3. 
The enhanced sensitivity available from a quadratic Hamil- 
tonian is a consequence of this increased rotation rate, which 
is greater by a factor of 2 J cos (3 — 2{J Z ) than that avail- 
able from a linear Hamiltonian. This same conclusion came 
out of the short-time analysis of Sec. IV, but it is stronger 
now because the uniform-fringe approximation is much better 
than the short-time approximation. The short-time approx- 
imation requires that J(f> <C 1 and thus describes correctly 
only the center of the central fringe for J y measurements. In 
contrast, the uniform-fringe approximation only requires that 
4>V~J\ sin/3| <C 1; within this requirement, there can be sev- 
eral fringes, i.e., 2 J <f) cos (3 can be somewhat larger than ir, 
provided that J 3> tan 2 (3. The more accurate uniform-fringe 
approximation allows us to see the other near-optimal oper- 
ating points for J y measurements and to see the optimal op- 
erating points for J x measurements, which lie not at <f> = 0, 
but at <ft = ±7r/4Jcos/3. As [3 approaches 7r/2, the fringes 
become wider and wider, making the uniform-fringe approx- 
imation reliable only for larger and larger values of J. For 
(3 = it/ 2, the fringes disappear entirely, and a separate analy- 
sis is required to find the scaling for J x measurements (since 
(Jy)^ — for [3 = it/2, measurements of J y provide no in- 
formation about <ft). 

That the final state (5.4), within the region of high sensi- 
tivity, is approximately an angular-momentum coherent state 
tells us two important things. First, even though the quadratic 
Hamiltonian will generate entanglement from a product state, 
this entanglement plays no role in the enhanced sensitivity. 
The improved sensitivity comes from the increased rotation 
rate of the coherent state, which is a product state, having no 
entanglement among the probe constituents. Indeed, for the 
measurements we consider here, the deviation from being a 
coherent state makes the sensitivity worse. Second, that the 
probe state is approximately a product state within the region 
of high sensitivity hints that this scheme should not be as frag- 
ile in the presence of decoherence as schemes that rely on ini- 
tial entanglement. We investigate the impact of decoherence 
in Sec. VB and show that the (D(n~ 3 / 2 ) scaling is unaffected 
by phase decoherence. 
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FIG. 1: (Color online) Sensitivity (solid red lines) vs. <j> (— tv/8 < 4> < n/8) using an optimal initial state at angle /3 = 7r/4: (a) J x 
measurements; (b) J y measurements. The total angular momentum J of the probe is 200, corresponding to n — 400. The lower bound on 
the sensitivity, 1/^/2^^ 2 , is plotted as the dotted (green) line. The sensitivity is characterized by rapidly oscillating fringes and a decay of 
sensitivity away from the best sensitivities near cj> — 0. The sensitivity patterns repeat with periodicity n; only a quarter of a period is plotted 
because the sensitivity worsens even more outside the plotted region. Part (a) also shows the sensitivity for J x measurements when — n/2 
(dashed blue line); notice the absence of fringes in this case and the substantially degraded sensitivity. 




FIG. 2: (Color online) Central few fringes of the measurement precision for /3 = n/4 and J — 2 500 (n = 5 000): (a) J x mea- 
surements; (b) J y measurements. The solid (red) lines are the exact sensitivities, the dashed (blue) lines are the sensitivities given 
by the uniform-fringe approximation of Eqs. (5.22) and (5.23), 

&<t>x ~ 2J3/2 ^ + sin 2 (2J^,co^y ' 54>v ~ 2J^{ 1 + cos 2 (2J^coi^y ' 
and the dotted (green) lines are the Gaussian-envelope approximation of Eqs. (5.17). The uniform-fringe approximation locates the fringes 
precisely, but misses entirely the degradation in sensitivity as one moves away from the central fringes and also fails to characterize accurately 
the shape of the fringes. The Gaussian-envelope approximation improves on this performance by capturing the degradation of sensitivity quite 
well, but still fails on the fringe shapes. Even the central fringe for J y measurements is noticeably flatter than in the two approximations. To 
get the best sensitivity, one should operate right on the central fringe, at <f> — qn, for J y measurements and on one of the two central fringes, 
centered at (j> — qn ± 7r/4Jcos/3 for J x measurements. Notice that J x measurements achieve nearly optimal sensitivity at points near the 
outside of these two central fringes. 



To achieve the optimal sensitivity for J x or J y measure- 
ments, we need to operate within the appropriate central 
fringe, of width 7r/v2J for {3 = 7r/4. This can be done by 
using an adaptive feedback procedure, which we discuss in 
the context of J y measurements. The feedback procedure is 
carried out in several steps, in each of which the quantity that 
is estimated is <fi — <f> cs t, where est is the estimate of <fi from 
the previous step. At each step, we choose J = n/2 so that 
4> — 4>cst is with very high probability close to the center of 
the central fringe, and we use v probes to determine cf> with 



greater precision for the next step. As we obtain progressively 
refined estimates of <j), the quantity being estimated becomes 
smaller and smaller, always lying well within a sequence of 
progressively finer central fringes. 

To check that this procedure works and to determine its 
scaling properties, imagine that we determine c^jl-n bit by bit. 
At step I, we determine the Zth bit of </)/2tt by choosing J = Ji 
so that the precision is given by 

84)i _ 1 1 1 
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where the factor / ~ 3-10 is chosen to ensure that we get the 
right kh bit with very high probability. This gives 



J i = 



1 



2i/V3 



2/3 



(5.28) 



We must, of course, choose J; to be an integer or half-integer, 
so we choose the nearest one, but this detail does not change 
the resource calculation significantly, so we ignore it. At 
step I + 1, 4> — 0cst lies well within the central fringe, as we 
see from 



f+i 



2 7/6 



1/3 



(5.29) 



Indeed, because of the 0(J~ 3 / 2 ) scaling, the quantity being 
estimated is buried progressively deeper fractionally in the 
central fringe as we step through the procedure, despite the 
fact that the central fringe is itself narrowing exponentially. 

Suppose now that we use this procedure to estimate L bits 
of <fi. The total number of constituents used, 



N = vJ2 2J i = 
i=i 



2vf\ 2/z 2 2L / 3 



1 



2 2 / 3 - 1 



(5.30) 



is dominated by the last step, as is typical in these feedback 
procedures. The ultimate precision displays the 0(iV~ 3 / 2 ) 
scaling, 



2tt2~ L = 



4/ 



(2 2 / 3 - If I 2 TV 3 / 2 

2/ 1 2 

(22/3 - 1)3/2 ^7 (JV/l/)3/2 



(5.31) 



with a small additional overhead given by the factor 
2//(2 2 / 3 - l) 3 / 2 . 

As j3 becomes smaller, the uniform-fringe approxima- 
tion becomes progressively better, since the fringes oscillate 
rapidly and the Gaussian envelopes become very broad. On 
the other hand, the signal in J x and J y disappears, making the 
sensitivity worsen as 1/ sin 2 2(3. 

At the other extreme, as (3 approaches tt/2, the uniform- 
fringe approximation becomes poorer as the fringes become 
as wide as the Gaussian envelopes and loses validity entirely 
when J | cot /?| ~ 1. When (3 = n/2, which is the initial probe 
state analyzed by Rey et al. [15] (Ref. [15] actually uses (3 — 
— tt/2, but this state is equivalent to (3 = tt/2 for purposes of 
these measurements), \J y )^ — 0, making J y measurements 
useless for extracting information about <f>. Thus we have to 
choose the J x measurement. The dashed line in Fig. 1 shows 
that the optimal operating point is cj> = (or, more generally, 
4> = qif). Near = 0, the expectation value and variance of 
J x are given by 



J(2J- 1) 

2 ' 

: J(2J-1)< 



(5.32) 
(5.33) 



where the approximations hold for i 
optimal sensitivity, 



< 1/ \fj . The resulting 



|d(J X ) /#| y/j(2J~l) 



(5.34) 



has the C(J _1 ) sensitivity scaling found in [15]. 

We can relate these results to the general lower-bound anal- 
ysis in Sec. Ill by noting that (3 = tt/2 means that (h) = 
(Z) /2 = 0. This means that the dominant sums in the ex- 
pansions of Eqs. (3.4) and (3.6) are those that contain only 
squares of hjS. The number of terms in these sums scales as 
0(n k ), which yields a sensitivity that scales as 0{n~ k / 2 ). 

To gain further insight into the scaling behavior, we plot 
the scaling exponent £ in 6(f) — 0(n~^) as a function of (3 
for J x measurements (Fig. 3) and J y measurements (Fig. 4), 
using three very large values of J. For J y measurements we 
calculate £ at the optimal operating point, = 0. For J x mea- 
surements, the optimal operating point is a function of j3, but 
a good compromise point, which works well over the entire 
range of/?, is 1/V2J, so we calculate the scaling exponent at 
this point for all values of /?. An investigation of nearby op- 
erating points scaling as 1/ J gives plots with no discernible 
differences for the large values of J under consideration. The 
main differences between J x and J y measurements are the 
following: (i) right at (3 = tt/2, J x measurements have a scal- 
ing exponent of 1 , whereas J y measurements provide no infor- 
mation about 7; (ii) for J y measurements, the plot of scaling 
exponent has two humps, nearly symmetric about (3 = tt/4 
and (3 = Ztt/A, whereas for J x measurements, the scaling 
exponent is better on the outside of the humps. The overall 
trend is for both measurements to have a scaling exponent of 
£ = 3/2 in the limit of large J, except at /? = 0, tt/2, and tt. 




FIG. 3: (Color online) Scaling exponent £ for J x measurements. The 
dotted (red) line is for J = 10 3 , the dashed (green) line for J = 10 5 , 
and the solid (blue) line for J = 10 7 . 



B. Decoherence 

The coherent-state model suggests that our generalized 
quantum metrology scheme with initial product states should 
not display the fragility of entangled protocols in the presence 
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FIG. 4: (Color online) Scaling exponent £ for J y measurements. The 
dotted (red) line is for J = 10^, the dashed (green) line for J = 10 5 , 
and the solid (blue) line for J — 10 7 . 



of decoherence. We can investigate this possibility by consid- 
ering independent dephasing of the effective qubits, described 
by the Lindblad equation 



-(ZpZ - p) , 



(5.35) 



where t 2 = is the dephasing time. Since dephasing com- 
mutes with the quadratic Hamiltonian, we can shunt its effects 
to the final time f, whence it maps the Pauli operators of each 
effective qubit in the following way: 



X 

Y -> e 
Z -> Z 



e- n X 
-rt. 



(5.36a) 
(5.36b) 
(5.36c) 



To obtain the effect of the decoherence on the expectation val- 
ues and variances of the measured operators at the time of 
measurement, it is easiest to use the adjoint map [3], which 
for this simple case is identical to the map (5.36) and gives 



{Jx,y)y 

(AJ Xt y)l = < 



-rt 



(Jx 



x,y/o 

-2Ttf\ j \2 



J 



(AJ^ + ^l-e"^) 



(5.37) 
(5.38) 



Here a subscript T denotes the value with dephasing, and a 
subscript without. It is now easy to see that under this model 
of decoherence, for either of the measurements that we are 
considering, the sensitivity takes the form 



hi 



<5 7 2 



S<f 1 + 



J(e 2Ft - 1) 
2v(d(J Xty ) /djy 
J(e 2rt -iy 



2(A Jx,y)o 



(5.39) 



To assess the effects of decoherence, we now focus on J, 



u 



measurements, and we assume that through an adaptive feed- 
back procedure like that sketched in Sec. V A, we are operat- 
ing well within the central fringe, i.e., jt is somewhat smaller 



than 7r/4Jcos (3. Inserting the = values from Eqs. (5.21) 
and (5.23) into Eq. (5.39) yields a sensitivity 



<57r 



£V^V2J 3 / 2 |sin2£| 



(5.40) 



If we now let T = vt be the total time available for measure- 
ments involving v probes, the optimal value of t, found by 



maximizing eF l j \ft is t = r 2 /2, gives a sensitivity 



^7r 



Tt 2 J 3 / 2 |sin2/?| 



(5.41) 



This result assumes that each probe can be processed in a 
time T2/2, but within this constraint, the scaling is the same 
0{ J~ 3 / 2 ) scaling that applies in the absence of decoherence. 
This is to be contrasted with entangled inputs, where uncor- 
rected phase decoherence degrades the scaling from 0(J~ 2 ) 
to the 0{ J~ 3 / 2 ) characteristic of product inputs. 

These arguments hold for general symmetric fc-body 
Hamiltonians, giving a sensitivity scaling 0(n~ k+1 ' 2 ) for ini- 
tial product states subjected to uncorrected phase decoher- 
ence. This is the same scaling achieved by initial optimal 
entangled states under this decoherence model [15]. On the 
other hand, the use of product states with fc-body Hamiltoni- 
ans for k > 2 can surpass both the standard quantum limit and 
the Heisenberg limit, even in the presence of phase decoher- 
ence. 



VI. CONCLUSION 

The possibility of using nonlinear Hamiltonians has the po- 
tential to open up a new frontier in quantum metrology. Quan- 
tum metrology has traditionally focused on linear Hamilto- 
nians of the form -fj z — T^Li^jA The main tech- 
nical challenge has been to improve on the standard quan- 
tum limit for determining the parameter 7, which scales as 
0(n" 1/2 ) and can be attained relatively easily using prod- 
uct input states and separable measurements. The goal of 
linear quantum metrology has been to achieve the Heisen- 
berg limit for determining 7, which scales as 0(7i _1 ) and 
requires the use of highly entangled input states. Nonlinear 
coupling Hamiltonians of the form offer the possibility of 
further improvements in scaling. With the same highly entan- 
gled input states, nonlinear Hamiltonians can achieve a scal- 
ing 0{n~ k ). More importantly, they provide 0(nr k+1 / 2 ) 
scalings, better than the Heisenberg limit, for input product 
states and separable measurements. We expect that the gen- 
eralized quantum metrology of nonlinear Hamiltonians will 
lead to new experiments — and, ultimately, to new devices — 
that take advantage of the enhanced scaling, which is available 
using the experimentally accessible tools of product-state in- 
puts and separable measurements. 

A notable feature of generalized quantum metrology is 
that the enhanced scalings available with product-state inputs 
do not rely on the entanglement produced by the nonlinear 
Hamiltonian. We reach this conclusion in this paper from a 
detailed analysis of the k — 2 case, in the course of which 
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we formulate an approximate coherent-state model of the time 
evolution, which applies during the period of enhanced sensi- 
tivity. In the model, a coherent state that makes an angle 
to the z axis rotates with angular velocity 2-fJ cos 0. The in- 
creased rotation rate, larger by a factor of 2 J cos than for 
k = 1, accounts for the enhanced sensitivity. Since coher- 
ent states are product states, this indicates that entanglement 
plays no role in the enhanced sensitivity, and it accounts for 
the robustness we find in the presence of phase decoherence. 

Although these conclusions emerge here from the k = 2 
analysis in this paper, it is not hard to extend the coherent- 
state model to arbitrary k. Given the input state (5.2), the state 
at time t = <f>/j becomes 

1*0(4)) - e-** I*,,) = ]T d m e^ mk | J, m) . (6.1) 

m 

The squares of the Wigner rotation-matrix elements d m of 
Eq. (5.3) are a binomial distribution, which for large J, ap- 
proaches a narrow Gaussian, centered at m = (J z ) = J cos 0, 
with half-width V2J + 1 sin 0. This encourages us to approx- 
imate m k in the phases of Eq. (6.1) as (J cos + Am) k ~ 
(Jcos0) k + k(J cos 0) k ~ 1 Am, giving 

\-qj p (t)} = e iMk-l)(Jco S f3) k J2 dme -m(Jco S f3)»-im jj^ 

m 

_ e i0O-l)(.7cos/3) fe e -i4>k{Jco S l3) k - 1 .] Ze -ip.J y | j 

(6.2) 

This is an angular-momentum coherent state at angle to 
the z axis, rotating about the z axis with angular velocity 
jk( J cos 0) k ~ 1 , which is the same enhanced rotation rate that 
we found in the very short-time analysis of Sec. IV. The ap- 
proximation leading to the coherent state (6.2) thus extends to 
arbitrary k the uniform-fringe approximation, formulated for 
k = 2 in Sec. V. The fringes have width ir/k(J cos 0) k ~ 1 , 
and the approximation provides a reasonable description of 
the first and second moments of J x and J y so long as 
(Jcos/3) fe - 2 ^>\/j| sin0\ < 1. 

The enhanced rotation rate is responsible for the improved 
scaling, and just as for k = 2, the coherent-state model indi- 
cates that the entanglement generated by the nonlinear Hamil- 
tonian plays no role in the enhancement. In separate work, to 
be published elsewhere, we extend these ideas. We investi- 
gate in more detail the entanglement generated by the nonlin- 
ear Hamiltonian, quantifying it using standard entanglement 
measures and showing that the enhanced sensitivity with ini- 
tial product states can be achieved with a vanishing amount of 
entanglement. 
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APPENDIX A: SYMMETRIC HAMILTONIAN WITHOUT 
SELF-INTERACTION TERMS 

The symmetric fc-body coupling without self-interactions is 
described by the Hamiltonian 

n 

H = ^ h a x ---h ak 
(ai,...,afe) 

n 

= k\ Y h ai --- K k 

a 1 <a 2 <-..<a k 



n 




where in the second form we use the expansions of Sec. Ill A 
to relate H to the Hamiltonian with self-interactions, plus cor- 
rections of which we give only the first. The analysis of the 
QCRB (1.6) relies on finding the largest and smallest eigen- 
values of H. Since we are only interested in the leading-order 
scaling of the QCRB, Eq. (Al) shows that the analysis pro- 
ceeds exactly as in the corresponding analysis for H in Sec. II, 
with the proviso that the results are only good to leading order 
in n. Thus for H, the QCRB (1.6), which can be achieved by 
using an appropriate initial entangled state, scales as 0(n~ k ). 

For the case of initial product states for the probe, we pro- 
ceed as in the comparable analysis of H in Sec. Ill A. Using 
the same trick of artfully switching between restricted and un- 
restricted sums, we can write 

(H)= Yl (h ai ) ■ ■ ■ (K k ) 

(ai,...,a k ) 

= Y M ■ ■ ■ (ha k ) 
ai,...,a k 

"(2) 51 { h a t ) ■ ■ ■ (K k _ 2 ) (h^^) 2 
(oi,...,a fc _i) 

+ 0(n k - 2 ), (A2) 



(H 2 )= £ (h ai )---(h a2k ) 

(a 1 ,...,a 2k ) 

(ai,...,02fe-i) 

+ 0(n k - 2 ), (A3) 
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(H) 2 = £ (h ai )---(h a2k ) 

(ai,...,a 2 k) 

+ k 2 ^ CO • • • ( h a 2k - 2 ) {K^f 

(«l,...,02t-l) 

+ 0{n k - 2 ) . (A4) 
The resulting variance of H, 

(AHf = k 2 ]T {K 1 )---{K 2k _,)Ahl k _ i 

(oi,...,a2js-i) 
2fc-2\ 



(A5) 



is the same, to leading order in n, as the variance (3.7) of H, 
so the remaining analysis of the optimal scaling and initial 
product states proceeds exactly as in Sec. III. 



APPENDIX B: FIRST AND SECOND MOMENTS OF J 

In this Appendix we derive the first and second moments of 
the components of the angular momentum J for the Hamilto- 
nian H y = 7 J| considered in Sec. V. 

Since J z commutes with the Hamiltonian, it is a constant of 
the motion, and its moments are at all times those of the initial 
coherent state: 

(J z ) 4> = Jcosp, (Bl) 
(j 2 )^ J 2 cos 2 /3+^sin 2 /3. (B2) 

In finding the first and second moments involving the equa- 
torial components of J, it is convenient to work in terms of 
the angular-momentum raising and lowering operators, J± = 
J x ±iJ y , which act according to J± \J,m) = T^j \ J, m ± 1), 
where T± = y/ (J T m)(J ± m + 1). 

For the evolved state (5.4), we can write 



m,m' — — J 

= £ d m d m+1 r+e^ 2m+1 >) 

771— — J 

J 

= cot(/3/2) ( J - ™)d 2 m e^ 2m+1) , (B3) 



m= — J 



where the last line uses 



d m +i — 



J — m 



- cot(/?/2K 



J + m + 

Using Eq. (5.3) and a derivative of the binomial formula, 

J 



(B4) 



£ (J-ro) 



m— — J 



2J 



J 



.,7-\-miJ—m 



2Jb(a + b) 



2J-1 



(B5) 



we obtain Eq. (5.9): 



<J + > = cot(/3/2) e * ]T V- m )[jl 

m— — J 

id) „„„2/ fl J+m r „~i<t> „;„2/a J - m 



2 J 

J — m 



x [e* cos 2 (/3/2)] J+m [e-* sin 2 (/3/2)] ' 
= Jsin/3(cos0 + isin(/)cos/3) 2j_1 . (B6) 

The evaluation of the remaining second moments proceeds 
along the same lines. The correlation between J z and J x or 
J y is conveniently expressed by 



- (J Z J+ + J+Jz)^ 

J 

= d m d m+1 (m+l/2)r+e^ 2m+i y) 

rn— — J 

= (1/2 - J) (J + ) 4> 

+ cot(/3/2) ]T (J 2 - m 2 )d 2 m e^ 2m +V 

m= — J 

— — — sin /3(cos (j) cos P + i sin </>) 



x (cos0 + zsin(/)cos/3) 2(j x) , 



(B7) 



where we use 



m— — J 



2J 



m~)\ \a J+m b 

J — m 



J+miJ—m 



2J(2J-l)ab(a + b) 



2(7-1) 



(B8) 



To find the second moments that involve only the equatorial 
components, we use 



J 2 + J 2 = -(J + J- + J-J + ), 

Jl-J 2 y=\(jl + J 2 -), 

J J + 7 / — — ( 7 2 - 7 2 ) 

<J x >->y \ f JyJx 2^ \ + — / ' 



(B9a) 
(B9b) 
(B9c) 



from which we get 
1 



(Jly), = j (J+J- + J-J+U ± ( J 2 ) ) , (BlOa) 



(J x Jy + JyJx)^ = Im( ( J 2 ) 4> ) . 

Thus we calculate Eq. (5.10), 



(BlOb) 



1 J 

-(J + J_+J_J + )^= Y, (J + J 2 -m 2 )d 



m= — J 

= J + 



J(2J ~ 1} sin 2 /? , (Bll) 
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where we use Eq. (B8), and we calculate Eq. (5.11), 

,7 

m— — J 

J 

= cot 2 (/3/2) ^ (J-m)(J-m- l)d 2 n e 4 ^ (m+1) 



m— — J 



J( - 2J - sin 2 /3(cos 2</> + i sin 2<j> cos /J) 2 ^ 1 ' 



(B12) 



where we use 



J2 (J-m)(J-m-l)( 

i— _ r \ 



m— — J 



2J 

J — m 



= 2J(2J-l)6 2 (a + 6) 2 < J - 1 ) 



(B13) 



The equatorial second moments listed in Eq. (5.16) and the 
cross moment, 



1 



1 



(J x J y + JyJJ^ = -Im( ( J+) ) 



= J(2J 4 1} sin 2 /? R 2 ( J -V sin[2(J - 1)6] , (B14) 

follow from inserting these results into Eqs. (BIO). 

We now make the uniform-fringe approximation of Sec. V, 
keeping only the fringe terms near = 0. This approximation 



requires that \fj <j>\ sin f3\ -C 1, but allows J<f> cos (3 to be con- 
siderably bigger than 1 when J cot 2 (3 is large. The resulting 
first and second moments, 



(J z ) = J cos/3 , 
(J+).~ Jsin/3 e 2 ^ c °^, 



(B15a) 
(B15b) 



< J 2 ) = J 2 cos 2 /? + J - sin 2 /? = J - + J(2J 2 1} cos 2 /3 , 

(B15c) 

i <J Z J+ + J + J z ) ~ ^-il sin/jcos/je 2 ^-^ , 

(B15d) 

- - 7 1 J ' 2J " 1) -sin 2 /3cos 2 (2J ( /)co S /3) , (B15e) 



a 2 ) - f + 
v'4 ~ 2 



( J , 2 ) <i ^ i + J{2 '\ 1) sin 2 ^sin 2 (2J0cos/3) , (B15f) 



'l^L — 11 sin 2 /? sin(2 J0 cos /?) cos(2 cos /?) 



(B15g) 

have the unique form of an angular-momentum coherent state. 
They show that in the uniform-fringe approximation, the state 
is an angular-momentum coherent state oriented at angle (3 to 
the z axis and rotating about the z axis with angular velocity 

27 J cos f3 = 2j(J z ). 
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